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Abstract 

We solve analytically the multiple scattering (KKR) equations for the two dimensional pho- 
tonic crystals in the long wavelength limit. Different approximations of the electric and magnetic 
susceptibilities are presented from a unified pseudopotential point of view. The nature of the 
so called plasmon-polariton bands are clarified. Its frequency as a function of the wire radius is 
discussed. 
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There is much interest recently in two dimensional (2D) photonic crystals (PC) consisting 
of arrays of metallic or dielectric cylinders (wires) in an insulating matrix or arrays of insu- 
lating cylinders in a metallic matrix. These include recent interest in left-handed materials[l] 
and in plasmonics. [5] A key issue is the effective susceptibilities < e > and < fi > of the 
system. To design systems at different frequencies such as in the infrared range it is useful 
to know their values for different system parameters. The photonic bands in an array of 
cylinders can be understood entirely in terms of the scattering phase shift of the cylinders. 
In the pseudopotential idea in electronic structure calculation a real potential is replaced by 
an effective one so that the same scattering of the electrons is produced. Similarly effective 
susceptibilities can be introduced so that the correct scattering effect for electromagnetic 
waves is produced. We examine this idea to derive effective susceptibilities for the cylinder. 
For example, from the scattering phase shift for a mode with the electric field along the 
cylinder axis, an effective dielectric constant for the cylinder is found to be 

e' E = -2J'e/(JkiR)[l + 0.5(k o R) 2 ln(k R)}/[l - fi kiRJ' ln(k R) / (J Hi)]; (1) 

e = 6j/e , J = Jo(kiR). The subscipts i, o refers to quantities inside and outside the cylinder, 
respectively. For metallic cylinders, when the skin depth is much less than the radius of 
the cylinder, the second term in the denominator is larger than the first term, we obtain an 
effective dielectric constant of a metallic form given by 

where the effective "plasma frequency" is given by 

co' p 2 = -2c 2 /[R 2 \n(coR/c)]; (2) 

R is the radius of the cylinder. While the original analysis [2] for the effective dielectric 
constant is carried out for a wire radius less than the skin depth, the experiments [3] for the 
left-handed materials are carried out for wires the width of which is larger than the skin 
depth. The above formula provides for an extension of the original analysis. There is a log 
correction term that has not been noticed before. 

This clarifies the issue of damping. For frequencies from 1 to 10 GHz, the imaginary 
part of the dielectric constant of most metals is about a thousand times larger than the real 
part. When the skin depth is much less than the wire radius, the loss in the metal is only 
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restricted near the surfaces of the wires and the effective damping is reduced. Indeed, the 
above effective dielectric constant depends only on the wavelength and the wire radius, with 
no damping! 

Eq (1) encompasses other recent results. For dielectric rods. If the second term of the 



denominator is smaller than the first term, we recover recent results by Wu et al. and Hu 
et al.fll, 12] that 



This pseudopotential idea is also implicit in recent results using cylinders with a high di- 
electric constant ferroelectric. [13] Eq. (1) extends these results to more general regions of 
the parameter space. 

In this paper we further calculate the photonic band structure of an array of cylinders 
of radius R in the long wavelength limit when the separation between the wires a is less 
than the free space wavelength A = 27r//c . We solve analytically the multiple scattering 
(KKR) equations in the long wavelength limit. We find that both the s and the p wave 
scattering phase shifts are of the same order of magnitude, (koR) 2 , and need to be considered. 
These produced two photonic branches, an "acoustic" mode with a frequency proportional 
to the wave vector with an effective dielectric constant < e > (eq. (1) and (8)) and a 
magnetic susceptibility < pi > (eq. (5) and (8)) and an "optic" mode with a gap. For 
negative susceptibilities and narrow cylinders, the "optic" mode corresponds to a flat band 
at frequencies close to the surface plasmon resonances, as has been previously discovered 
numerically. For the acoustic mode, we found that < e > can be expressed as the arithmetic 
mean of that of the medium and an effective dielctric constant of the cylinder e' c . We now 
describe our result in detail. 

Pseudopotential: We first describe our "pseudopotential" idea for the effective dielec- 
tric constant of the cylinder. As far as the EM field outside the cylinder is concerned, all 
that matters is the scattering phase shift for angular momentum component n given, for the 
E (TM) mode, by 



kj = ko(fijej) 1 / 2 for j=o,i, e = ej/e . Similarly, for the H (TE) mode ki = (ej/ij) 1 / 2 /^. 




2J'e/(JkiR) 




J'^koRjkiJnjkiR) - J' n {kjR)Jn{k R)k e 
kiJ n (kiR)N^k R) - k eN n {k R)J' n {hR) 



(3) 



tan 7], 



'n 



H 



J n {k R)k 6J n (kiR) J' n (kiR)) J n (k R)ki 
k eJ n {hR)N' n (k R) - kiN n {k R)J' n {kiR)' 



(4) 
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We first focus on the s wave with n=0. For k R « 1, J = 1, J' Q = —k Q R/2, N = 
2ln(k R)/ir, = 2/(Txk R). 

tan 77^ « -0.25(A; OJ R) 2 7r[l + 2J'e/{Jk i R)]/[l - ^hRJ' 'ln(k R) / '(J 

tan^ = -0.25(A; OJ R) 2 7r[l + 2J'ife i /( JeA; 2 i?)]/[l - k i R\n(k R)J l /(Je)}. 
For fcji? also small 

tan 77^ Rj — 7r(fc i?) 2 [l - e]/4; tan 77^ rj -ir(k R) 2 [l - /ij//i ]/4 

As is expected, when e = 1, there is no scattering and tan 77^ = 0. 

When kiR is not small, one can define effective susceptibilities so that the same phase 
shift is produced: 

tan77 £ « -7r(k R) 2 [l - 4]/4; tan 77^ « -7r(A; i?) 2 [l - /4//W 4 (5) 

This is the "pseudopotential" idea that we mentioned. From eq. (5) we obtain eq. (1) for 
the effective dielectric constant and also an effective magnetic susceptibility for the H mode: 

H' H = -[kiR\n(k R)J'/(Je) + 2J'k i /{Jek 2 R)]/[l - kiR\n(k R)J'/(Je)}. (6) 

The second term in the denominator is of the order of (fii/ei) ' 5 k o R and is usually smaller 
than the first term. We obtain 

/4 « -2J'ki/(Jtk 2 R). 

We next investigate the phase shifts for the higher order partial waves. In the limit 
k a R « 1, 

tan 77^ = ix(k R/2) 2n / \{n - l)!n!)[/i - ^nJ n /{ J^R)}/^ + ^nJJiJ^R)}, 

tanr7 T f = ir(k R/2) 2n /[(n - l)!n!][l - neJ n /(J' n kiR)]/[l + neJj (JfaR)]. 
Here J n = J n {kiR). When kiR is also small 

tan 77^ = [(k R/2) 2n /n}(fi-l)/(fi+l); tan 7^ = [(k R/2) 2n /n](e - l)/(e + 1). 

There is recently much interest in "plasmonics" when the frequency is close to the interface 
plasmon frequency so that e = — 1.. At this frequency 77^ = tx/2. A Mie scattering resonance 
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is exhibited for the TE modes for all n ^ 0. For n = 1 the requirement that the same 
scattering phase shift is obtained even when kiR is not small provides for the equations 
determining the effective susceptibilities: 

tariff = -n{k R/2) 2 {fi - fi' E )/{fi' E + fi Q ). (7) 

tan 77? = -n(k R/2) 2 (e - e' H )/(e + e' H ). (8) 
From these we obtain the effective susceptibilities 

/4 = HiJi/(J[kiR); e' H = eiJJiJikiR). (9) 

Similar equations have also been obtained by Hu et al. and Wu et al.fll, 12] from a coherent 
potential approximation. The results here provides a different interpretation of their results. 
With the current view, "plasmonics" phenomena can also be manifested for no n- metallic rods 
if e' H + e Q = 1 and the same scattering phase shift is produced. We next turn our attention 
to the photonic bands. 

Photonic band structure: First we briefly recapitulate the multiple scattering (KKR) 
technique. The basic idea is that the scattered wave from the photonic crystal is self- 
consistently sustained. More precisely, consider a cylinder at the origin. The scattered wave 
from all the other cylinders sum to produce a net incoming wave at the origin which is 
scattered by the cylinder at the origin and in turn produce a scattered wave from the origin. 
This scattered wave from the origin is related to the scattered wave from the other sites by 
a phase factor determined by the wave vector. 

More precisely, we denote the amplitude of the partial scattered wave with angular mo- 
mentum n by a n . The sum of the scattered waves from all the other sites becomes an 
incoming wave at the orgin with the amplitude p n = J2 n a n'S{n — n') where 

S(m) = ^ exp(ik ■ H)H m (k R) exp(im(f) R ). (10) 

Note that S does not include the wave from the origin; thus the sum in eq. (10) does 
not include the term at R=0. The outgoing scattered wave at the origin is related to the 
incoming wave by the t matrix: a n = t n p n . Substitute in the definition of p n , we arrive at 
the equation det[S(n — n') — 5{n — n')/t n ] = 0. Since the t matrix is related to the phase 
shift by t n = tanr] n /(tan?7 n + i), we obtain the KKR equation 

det[A(n - n') - S(n - n') cot rj n ] = 0. (11) 
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Here the structure factor A(m) = [S(m) — In this paper, we have assumed a time 
dependence of exp(iut). Outgoing spherical waves correspond to a H n = J n + iN n . 

In the long wavelength limit, one can approximate the sum for S by an integral, which 
can then be analytically evaluated. [14] The structure factor becomes 

A(n) » 4[ieMi<f>k)k/k ] n /[(k a) 2 (k 2 /k 2 - 1)] 

As is discussed above, if the wavelength outside the cylinder is long and k Q R << 1, tan^„ oc 
(k R) 2n for n^O, tan^ oc (k R) 2 . The s and the p wave shifts are of the same order of 
magnitude, (k R) 2 , and need to be considered. When only the s and the p wave components 
are kept, the KKR equation now becomes 



HE = 



(12) 



where 



H 



(13) 



A(0) + cot??i A(l) A(2) 
A(l)* A(0) + cotr] A(l) 
A(2)* A(l)* A(0)+cot?7i 

There are two classes of solutions, with either E x = E*_ x = \Ei\ exp(i(j) k ) or E 1 = —E*_ x 
i\Ei \ exp(i(f) k ). We get two possible eigenvalue equations. The first one is given by 

[A(0) + cot Vl + \A(2) \] [A(0) + cot r] ] - 2\A(1)\ 2 

For the second case, we get 

A(0) + cot77i + |A(2)| =0. (14) 

As we show below, the first mode corresponds to an "acoustic" branch with a frequency 
proportional to the wave vector with effective susceptibilities for the system; the second 
mode corresponds to a band with a gap. For negative susceptibilities, this corresponds to 
a flat band at frequencies close to the surface plasmon resonances, as has been previously 
discovered numerically [9]. 

We discuss the acoustic branch first. Substituting in the expressions for the phase shifts 
(eq. (5) and (7)) and the structure factor into eq. (11) and after some algebra, we obtain 



/,•- 



k 2 < e >< fi > 



where 



< e >= e G (l - /) + fe h 



(15) 



< /i >= M{\ +f) + fi (l - /)]/[/4(l - /) + /i (l + /)]. (16) 

In the static (zero wavevector and frequency) limit [4] for the case with the E field along the 
axis the tangential component of the electric field in the cylinder (i) and the medium outside 
(o) is the same: E Q = Ei = E. The average displacement field is given by < D >= c D + 
(1 — c )Di. Since D Q = e E , Di = e^Ei, we obtain < D >= (c e + (1 — c )ei)E =< e > E. 
Hence the average dielectric constant is just the arithmetic mean of the dielectric constants 
of the components: < e >= (c e + (1 — c )ej). This is the same as eq. (15). 

In multilayer systems, a similar result is obtained. [8] In that case the effective \x is the har- 
monic mean of the components while the effective dielectric constant is still the arithmetaic 
mean of that of its components. 

We next discuss the "optic" mode. Substituing in the expressions for the phase shifts 
and the A's, the equation for the second optic mode becomes 

2 In (k a/2^) (k a) 2 /n = 4 - 4/- 1 (/i c + /4)/(/4 - n ) + 0(k 2 ) 

for the E mode and 

2 In (k Q a/2^) (k af/n = 4 - 4r 1 (e D + e' H )/{e' H - e a ) + 0(k 2 ) 

for the H mode. When k approaches zero, k Q is not zero! Let us illustrate the physics by 
looking at the H mode. The limit of small f is particularly interesting. In that limit, the 
frequency is determined by the condition that e G + e' H = where e' H is given in eq. (8). 

For metallic cylinders with their radii less than the skin depths, hiR « 1, e' H — — 
1 — ujp/u 2 where u p is the plasma frequency. e Q + e« = when oj is equal to the inteface 
plasmon resonance, cu sp = u) p /(l + e^ 1 / 2 , For small k, from the above equation, we see that 
uj{k) = u){k = 0) + O(fk). If / is small, the dispersion is weak and the bands are flat. This 
flat band has been observed numerically previously. [9] The present calcualtion provides for 
a more direct analytic demonstration of this result. If k{R is not small, eq. (8) suggests that 
even with insulating cylinders, flat "plasmonic" photonic bands can still be obtained if the 
condition 

tiJiKJ'^iR) = -e . 

is satisfied. 

Let us next look at the E mode, the condition becomes // E /// = — 1- We call this the 
"magnetic surface plasmon" mode. Even though there is a lot of interest in plasmonics that 
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focus on the condition e' H /e Q = — 1, the corresponding condition on fj, have not been much 
discussed. 

There is another way to think of this type of solutions. As can be seen from eq. (5) 
and (6), when the susceptibilities of the metal are negative, e Q + e m can become zero and 
tan 7] = oo. The scattering can go though resonances due to the interface plasmon. This 
can lead to flat photonic bands, as has been observed in previous numerical calculations. In 
general, the more rapidly varying the phase shift, the flatter the band. 

Pokrovsky and Efros[10] have recently investigated the propagation of electromagnetic 
(EM) waves in a periodic array of metallic cylinders (wires) in the limit kR » 1. Our 
conclusion differs from theirs. In their paper, an expression similar to our 5(0) also appears. 
However, in their expression, the sum is over all R whereas the R=0 term is absent in ours. 
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The dominant contribution in the long wavelength limit is obtained by replacing H by iN. 
A(n) « 4[iexp(i<f> k )k/k ] n [l + 0{k 2 a 2 )]/[(k a) 2 {k 2 /k 2 - 1)] 
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